#ifndef INVERT_H
#define INVERT_H

#include <math.h>
#include <stdio.h>

void invert_matrix ( const double A[4][4], double Ainv[4][4] );

void adjoint( const double in[4][4], double out[4][4] );
double det4x4( const double m[4][4] );

double det3x3( double a1, double a2, double a3, double b1, double b2, double b3, double c1, 
	double c2, double c3 );

double det2x2( double a, double b, double c, double d );


#endif

